Near-term forecasting of Covid-19 cases and hospitalisations in Aotearoa New Zealand (preprint)

Near-term forecasting of infectious disease incidence and consequent demand for acute healthcare services can support capacity planning and public health responses. Despite well-developed scenario modelling to support the Covid-19 response, Aotearoa New Zealand lacks advanced infectious disease forecasting capacity. We develop a model using Aotearoa New Zealand's unique Covid-19 data streams to predict reported Covid-19 cases, hospital admissions and hospital occupancy. The method combines a semi-mechanistic model for disease transmission to predict cases with Gaussian process regression models to predict the fraction of reported cases that will require hospital treatment. We evaluate forecast performance against out-of-sample data over the period from 2 October 2022 to 23 July 2023. Our results show that forecast performance is reasonably good over a 1-3 week time horizon, although generally deteriorates as the time horizon is lengthened. The model has been operationalised to provide weekly national and regional forecasts in real-time. This study is an important step towards development of more sophisticated situational awareness and infectious disease forecasting tools in Aotearoa New Zealand.

Improving estimates of epidemiological quantities by combining reported cases with wastewater data: a statistical framework with applications to COVID-19 in Aotearoa New Zealand (preprint)

Background: Timely and informed public health responses to infectious diseases such as COVID-19 necessitate reliable information about infection dynamics. The case ascertainment rate (CAR), the proportion of infections that are reported as cases, is typically much less than one and varies with testing practices and behaviours, making reported cases unreliable as the sole source of data. The concentration of viral RNA in wastewater samples provides an alternate measure of infection prevalence that is not affected by human behaviours. Here, we investigated how these two data sources can be combined to inform estimates of the instantaneous reproduction number, R, and track changes in the CAR over time. Methods: We constructed a state-space model that we solved using sequential Monte Carlo methods. The observed data are the levels of SARS-CoV-2 in wastewater and reported case incidence. The hidden states that we estimate are R and CAR. Model parameters are estimated using the particle marginal Metropolis Hastings algorithm. Findings: We analysed data from 1 January 2022 to 31 March 2023 from Aotearoa New Zealand. Our model estimates that R peaked at 2.76 (95% CrI 2.20, 3.83) around 18 February 2022 and the CAR peaked around 12 March 2022. Accounting for reduced CAR, we estimate that New Zealand's second Omicron wave in July 2022 was similar in size to the first, despite fewer reported cases. We estimate that the CAR in the BA.5 Omicron wave in July 2022 was approximately 50% lower than in the BA.1/BA.2 Omicron wave in March 2022. The CAR in subsequent waves around November 2022 and April 2023 was estimated to be comparable to that in the second Omicron wave. Interpretation: This work on wastewater-based epidemiology (WBE) can be used to give insight into key epidemiological quantities. Estimating R, CAR, and cumulative number of infections provides useful information for planning public health responses and understanding the state of immunity in the population. This model is a useful disease surveillance tool, improving situational awareness of infectious disease dynamics in real-time, which may be increasingly useful as intensive pandemic surveillance programmes are wound down.

Simulating the impact of vaccination rates on the initial stages of a COVID-19 outbreak in New Zealand (Aotearoa) with a stochastic model (preprint)

Aim: The August 2021 COVID-19 outbreak in Auckland has caused the New Zealand government to transition from an elimination strategy to suppression, which relies heavily on high vaccination rates in the population. As restrictions are eased and as COVID-19 leaks through the Auckland boundary, there is a need to understand how different levels of vaccination will impact the initial stages of COVID-19 outbreaks that are seeded around the country. Method: A stochastic branching process model is used to simulate the initial spread of a COVID-19 outbreak for different vaccination rates. Results: High vaccination rates are effective at minimizing the number of infections and hospitalizations. Increasing vaccination rates from 20% (approximate value at the start of the August 2021 outbreak) to 80% (approximate proposed target) of the total population can reduce the median number of infections that occur within the first four weeks of an outbreak from 1011 to 14 (25th and 75th quantiles of 545-1602 and 2-32 for V=20% and V=80%, respectively). As the vaccination rate increases, the number of breakthrough infections (infections in fully vaccinated individuals) and hospitalizations of vaccinated individuals increases. Unvaccinated individuals, however, are 3.3x more likely to be infected with COVID-19 and 25x more likely to be hospitalized. Conclusion: This work demonstrates the importance of vaccination in protecting individuals from COVID-19, preventing high caseloads, and minimizing the number of hospitalizations and hence limiting the pressure on the healthcare system.